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\ abstract 

We prove factorization theorem for the processes 7r7* — > 7 and itj* — ► 7r to leading twist in the covariant gauge by 
I means of the Ward identity. Soft divergences cancel and collinear divergences arc grouped into a pion wave function 
£\j ■ defined by a nonlocal matrix element. The gauge invariance and universality of the pion wave function are confirmed. 
^ ' The proof is then extended to the exclusive B meson decays B — ► 7ZP and B — > 7rZi? in the heavy quark limit. It is 
shown that a light-cone 2? meson wave function, though absorbing soft dynamics, can be defined in an appropriate 
frame. Factorization of the B — > irlv decay in kx space, kx being parton transverse momenta, is briefly discussed. 
We comment on the extraction of the leading-twist pion wave function from experimental data. 
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I. INTRODUCTION 



The fundamental concept of perturbative QCD (PQCD) is factorization theorem, which states that nonperturbative 
dynamics of a high-energy QCD process either cancel or can be absorbed into hadron wave functions. The remaining 
part, being infrared finite, is calculable in perturbation theory. A full amplitude is then expressed as the convolution of 
a hard amplitude with hadron wave functions. A wave function, because of its nonperturbative origin, is not calculable. 
However, PQCD still possesses a predictive power, since a wave function is universal, i.e., process-independent. With 
this universality, a wave function determined by some means, such as QCD sum rules and lattice theory, or extracted 
from experimental data, can be employed to make predictions for other processes involving the same hadron. 

Nonperturbative dynamics is reflected by infrared divergences in radiative corrections. There are two types of 
infrared divergences, soft and collinear. Soft divergences come from the region of a loop momentum I, where all its 
components vanish: 

i** = (i+ r,Zr)~(A,A,A) . (1) 

The light-cone variables have been adopted, and A is a small scale. Collinear divergences are associated with a massless 
quark of momentum P = (Q, 0, 0t), Q being a large scale. In the collienar region with I parallel to P, the components 
of I behave like 

l^(Q,X 2 /Q,\) . (2) 

In both regions the invariant mass of the radiated gluon diminishes as A 2 , and the corresponding loop integrand may 
diverge as 1/A 4 . As the phase space for loop integration vanishes like d^l ~ A 4 , logarithmic divergences are generated. 

Factorization of the above infrared divergences in a QCD process needs to be performed in momentum, spin, and 
color spaces. Factorization in momentum space means that a hard amplitude does not depend on the loop momentum 
of a soft or collinear gluon, which has been absorbed into a hadron wave function. Factorization in spin and color spaces 
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means that there are separate fermion and color flows between a hard amplitude and a wave function, respectively. 
To achieve these, we rely on the eikonal approximation for loop integrals in leading infrared regions, the insertion of 
the Fierz identity to separate fermion flows, and the Ward identity to sum up diagrams with different color structures. 
Under the eikonal approximation, a soft or collinear gluon is detached from the lines in a hard amplitude and in other 
wave functions. The Fierz identity decomposes a full amplitude into contributions characterized by different twists. 
The Ward identity is essential for proving factorization theorem in a nonabelian gauge theory. 

In this paper we shall derive the factorization of the scattering processes 7T7* — * 7 and 717* — * tt, which involve the 
pion transition form factor and the pion form factor, respectively, at leading twist using the above techniques. Infrared 
divergences, occuring at higher powers of 1/Q 2 , Q being momentum transfer in the above processes, are neglected. 
It will be shown that soft divergences cancel and collinear divergences, factored out of the whole proesses order by 
order, are absorbed into a pion wave function, which is defined by a nonlocal matrix element. The universality of the 
pion wave function is equivalent to the universality of the collinear divergences in the two processes. The definition of 
the hard amplitude at each order will be given as a result of the proof. We emphasize that our derivation is simple, 
explicitly gauge- invariant, and appropriate for both the factorizations of the soft and collinear divergences. 

The factorization of the process 717* — > has been proved in jjj, but in the axial (light-cone) gauge A + = 0. 
In this gauge factorization automatically holds and the analysis is straightforward, because collinear divergences 
exist only in reducible diagrams. Our proof is performed in the covariant gauge, in which collinear divergences also 
exist in irreducible diagrams. The colleciton of these collinear gluons forms a path-ordered integral along the light 
cone, which renders a hadron wave function explicitly gauge invariant. The pion wave function has been constructed 
from 7*7 — > tt in the framework of covariant operator product expansion However, it was not clear how the 

path-ordered integral is generated We shall demonstrate that it appears as a consequence of the Ward identity. 
The factorization of 7T7* — > 7r has also been proved in Q based on a generalization of the Zimmermann's "A- forest" 
prescription which involves complicated diagram subtractions. It will be found that our derivation is simpler. 

We then prove factorization theorem for the exclusive B meson decays B — > jlu and B — > ttIv, whose topologies 
are similar to the scattering processes 7T7* — > 7 and 7T7* — > tt, respectively. The infrared divergences in B — > ^Iv 
and B — » ttIv have been investigated in || and [Q, respectively. However, a rigorous proof of factorization theorem 
has not yet been available. In the heavy quark limit, terms of 0(Aqct>/Mb), Aqcd being the QCD scale and Mb 
the B meson mass, are higher-twist and negligible. It will be shown that in this limit a gauge-invariant B meson 
wave function and hard amplitudes can be defined. The B meson wave function absorbs soft divergences of the above 
decays, which differ from the collinear divergences in the pion wave function. However, it is still possible to construct 
a light-cone B meosn wave function Q , if an appropriate frame with the photon (pion) moving in the minus or plus 
direction is chosen. 

We investigate the 0(a s ) collinear divergences contained in the process 7T7* — > 7 in Sec. II, and present the all-order 
proof of its factorization theorem by means of the Ward identity in Sec. III. The factorization of the process 7T7* — > tt 
is derived in Sec. IV. In Sec. V we prove factorization theorem for the B — ► 7ZP and B — > ttIv decays. Section VI 
is the conclusion. The detailed evaluation of the 0{a s ) collinear divergences in 7T7* — * tt are supplied in Appendix 
A. In Appendix B we comment on the extraction of the leading-twsit pion wave functions from experimental data of 
7T7* — > 7 and 7T7* — > tt. 

II. O(as) FACTORIZATION OF ttj* -> 7 

We start with the factorization of the process 7T7* — > 7 at the one-loop level, which will serve as the basis of the all- 
order proof. The momentum Pi of the pion and the momentum Pi of the outgoing on-shell photon are parametrized 
as 

P 1 = (P+,0,O T ) = -^(l,0,O T ), 

P 2 = (0,P 2 -,O r ) = -^=(0,l,0 T ). (3) 

Let e denote the polarization vector of the outgoing photon, which contains only the transverse components. We 
consider the kinematic region with large Q 2 = —q 2 , q = P^ — P\ being the virtual photon momentum, where PQCD 
is applicable. 

The lowest-order diagrams are displayed in Fig. 1. The lower valence quark (an anti-quark) in the pion carries 
the fractional momentum xP\. Contracting the amplitude in Fig. 1(a) with the leading-twist structure /Pi 75/ \/2Nc 
associated with the pion, N c — 3 being the number of colors, we derive 
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r (o),^ _ „-„2. jN c tr[jti{jP 2 - x A)7m A75] _ ,„ 2 . /^c ir(/ A75) 



lj V2 (P2 — a; Pi) 2 V2 xQ 2 ' () 

Figure 1(b) leads to the same amplitude, because the pion wave function is symmetric under the interchange of x and 
1 — x. The internal quarks are regarded as being hard, i.e., being off-shell by 0{Q 2 ), since contributions from the 
small x region will be suppressed by the pion wave function, which vanishes like x as x — ► 0. 

We identify the infrared divergences from 0(a s ) radiative corrections to Fig. 1(a), which are shown in Fig. 2. 



Self-energy correction to the internal quark, giving a next-to-leading-oreder hard amplitude, is not included. Figures 
2(a)-2(c) are the reducible diagrams with the additional gluon attaching the two valence quarks of the pion. It has 
been known that soft divergences cancel among these diagrams. The reason for this cancellation is that soft gluons, 
being huge in space-time, do not resolve the color structure of the pion. Collinear divergences in Figs. 2(a)-2(c) do 
not cancel, since the loop momentum flows into the internal quark line in Fig. 2(b), but not in Figs. 2(a) and 2(c). 
To absorb the collinear divergences, we introduce a nonperturbative pion wave function. 

The factorization of Figs. 2(a)-2(c) is achieved by meas of the insertion of the Fierz identity, 

hjhk = \hkhj + j(75)*(75)ij + ^h a )tk(l a )ij + ^(757ak-(7 Q 75)ij + ^{° a f})ik{o a(S )ij , (5) 

where / represents the identity matrix, and a a p is defined by a a p = i["f a ,^//3]/2. Different terms in the above identity 
lead to contributions of different twists. Take Fig. 2(b) as an example, whose loop integrand is given by 

r 2 2 r [N~c, \ x /Pi- / i fo-x A+ / 
h b = e g C F ^-tr^ V{xPi _ l)2 ^—^—^ 

(1 - x) /Pi+ / „ \ 1 

X [(l-x)P 1+ l}^ ^}p> (6) 

with Cf being a color factor. Inserting the Fierz identity, we obtain 

x Pi- / 757a (1 ~ x) /Pl+ / „ 7~75 ] 1 
{xPx-l) 2 2 [(l-x)Pi-M] 27 2 [P 



hb ~ ig C F tr{ 7„ 



/ ■ 2 ^ [K tr[/t{fo - x Pi+ /) 7At 7"75)Pi + (7 , 

xHe, VT (P2-XP1 + 2 ■ (7) 

Obviously, in the collinear region with the loop momentum I parallel to Pi, only the pseudo- vector structure 757a 
contributes to the first trace as shown above. Moreover, the matrices j a and 7" must be 7_ = 7 + and 7", respectively. 
Equation ([?]), as integrated over /, is rewritten as the convolution of the lowest-order hard amplitude H^(^) with 

the 0(a s ) pion wave function <p^ in the momentum fraction £ = x — l + / P^: 

/ al «^ (0) (0. (8) 

CFtr \^Jxl^lf^ [(l-x)P 1+ l}^ —}P ■ (9) 

(p2b contains the collinear divergence in Fig. 2(b), because the integrand in Eq. (0) diverges as 1/A 4 . The dependences 
on l~~ and on It in H^°\ being subleading according to Eq. (||), have been neglected. 

Diagrams with the additional gluon attaching the internal quark, Figs. 2(d) and 2(e), do not contain soft divergences, 
because the internal quark is off-shell. For example, the loop integrand corresponding to Fig. 2(d) is approximated, 
in the I — ► region, by 

' ' 0(A- 3 ), (10) 
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(p 2 - xP x + o 2 [(i - x)p l + z] 2 / 2 2(1 - x){p 2 - X p l yp l ■ u 

which is suppressed by the phase space for loop integration d I ~ A 4 . Therefore, we consider only the collinear 
divergences. The integrand associated with Fig. 2(d) is written as, 
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hd = -e z g 2 C F \ Hrtrl £ 



Pi-x p x Pi-x p x + / 



{p 2 - x p 1 y lv (p 2 - x p 1 + i) 



-If 



1 

72 



(11) 



Since / and 7 M involve only 7^, we drop —x /P\ and — x /Px + / in the internal quark propagators in the collinear 
region. Equation (0) is then simplified into 
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e 2 g 2 C F 



2P 2i 



2 (P2-ZP1 + Z) 



2 1 ^(P 2 ~xP 1 ) 27 ' 1 [(l-a ; )P 1 +;] 27 ^ l7 ° 



For the /-dependent hard propagator, we employ the relation 



2P 



2// 



(P 2 - .xPi + If 



n ■ I 



(P2-ZP) 2 



(p 2 - xPi + o 2 



(12) 



(13) 



which is an example of the Ward identity. The dimensionless vector n = P 2 /P 5 7, representing the direction of an 
cikonal line, lies along the outgoing photon momentum. n v is called the eikonal vertex and 1/n ■ I is called the eikonal 
propagator. The appearence of the eikonal line is a consequence of the Ward identity. Inserting the Fierz identity, we 
derive the factorization, 



with the 0(a s ) pion wave function 
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-ig C F tr< 



757+ (1 - x) Pi+ / 7 75 1 1 n v 



2 [(1 - x)P x + 1} 



r7 



Pn-l 



(14) 



(15) 



The collinear contribution from Fig. 2(d) has been split into two terms, with the first and second terms described by 
Figs. 3(a) and 3(b), respectively, where the eikonal propagator is represented by a double line. In Fig. 3(b) the loop 
momentum I flows into the internal quark line, such that the second term is a convolution of with (f)^) in the 
momentum fraction £. 

Figure 2(e) gives the loop integrand 



he = e'g'C F \ Mtr 7 



v x Px- / .P2-XP1+/ 



xPx-if ^(p 2 -xP 1 + iy 



P2-X p x 



If A 75 



(P 2 -xP 1 f 

Following the similar procedures, we obtain, in the collinear region, 

7 2e «-4eW°^)-ff (0) (*)], 

with the 0(a s ) pion wave function 

v A- / 7s7 + 7~75 1 1 n v 



ig z C F tr{ 7" 



(xPx - If 



Pn-l 



(16) 



(17) 



(18) 



The first and second terms in Eq. (|i7| ) are described by Figs. 3(c) and 3(d), respectively. 

Comparing Eqs. (|^), ( |l5| ) and (|18j), the Feynman rules for the perturbative evaluation of the pion wave function 
are clear: eff 1 ^ can be written as a nonlocal hadronic matrix element with the structure 7s7 + /2 sandwiched, which 
comes from the insertion of the Fierz identity: 



(j) {1 \x) 



/2lVcP 1 + 



(19) 
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The sum over colors is understood. The integral over z in fact contains two pieces: For the upper eikonal line in 
Fig. 3(a), z runs from to oo. For the lower eikonal line in Fig. 3(b), z runs from oo back to y" . It is easy to confirm 
that the above definition reproduces all the leading-twist collinear divergences in Figs. 2(a)-(c) and in Figs. 3(a)-3(d). 
The light-cone coordinate y~ ^ corresponds to the fact that the collinear divergences in Fig. 2 do not cancel exactly. 
At last, the factorization formula for the process nj* — > 7 is written, up to 0(a s ), as 

(0(0) + 0(1)) g, (i3 -(0) + + Q( a 2} ^ (2Q) 

with </>(°' = 1 and <g) representing the convolution in the momentum fraction. The 0(a s ) hard amplitude H^- 1 ' is 
defined by, 

tf (1) -E/(^-^ (1) ^ (0) > ( 21 ) 

where the subscript i runs from 2a to 2e. Obviously, is infrared finite, since all the (a s ) collinear divergences 
have been absorbed into the pion wave function (f>^ . 



III. ALL-ORDER PROOF OF FACTORIZATION THEOREM 



In this section we present the all-order proof of leading-twist factorization theorem for the process 7r7* 
construct a gauge-invariant pion wave function defined by 



1 



dy 
2ir 



e^ + MO|g(y-)^_Pexp 



-ig I dzn ■ A{zn) 
'0 



k(Pi)>, 



7, and 



(22) 



as shown in Fig. 3(e). The notation P means the path ordering. By expanding the quark field q(y~) and the path- 
ordered exponential into powers of y~ , the above matrix element can be expressed as a series of covariant derivatives 
(D + ) n q(0) , implying that Eq. ( p2|) is gauge invariant. 

It has been mentioned in the Introduction that factorization of a QCD process in momentum, spin, and color spaces 
requires summations of many diagrams, especially at higher orders. Hence, the diagram summation must be handled 
in an elegant way. For this purpose, we employ the Ward identity, 



1^(1, k u k 2 , 







(23) 



where G M represents a physical amplitude with an external gluon carrying the momentum I and with n external 
quarks carrying the momenta ki, k 2 , ■ ■ •, k n . All these external particles are on mass shell. The Ward identity can 
be easily derived by means of the Becchi-Rouet-Stora (BRS) transformation pl| . 

We prove factorization theorem by induction. The factorization of the 0(a s ) collinear divergences associated with 
the pion has been worked out in Sec. II. Assume that the factorization holds up to O(a^): 



G 



H 



with 



A' 



(24) 



(25) 



i=0 



i=0 



i=0 



G (l) denotes the full diagrams of 0{a\) with G (0) = 
0(a\) terms in the perturbative expansion of Eq. (t 
then have the relations 



H<® in Eq. (|), the pion wave function cj>^ is defined by the 
a), and flw is the infrared-finite hard amplitude of 0(a l s ). We 



G (i) = 



3=0 



(26) 



for i = 0, • • • , N, which imply that all collinear divergences in G^ have been collected into ^ , 3 ^ h systematically. 
Consider a complete set of full diagrams 

Gr(w+i) of o(a* +1 ). We look for the gluon, one of whose ends attaches the 
outer most vertex on the upper quark line in the pion. Such a gluon exists, since 

G (N+1) 

are finite-order diagrams. 
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We then classify G^ N+1 ^ into two categories, reducible and irreducible, according to the attachment of the other end 
of this gluon. If the other end attaches the upper or lower quark lines directly, the diagrams are reducible. The 0(a s ) 
examples are Figs. 2(a)-2(c). All other diagrams, with the other end attaching inside of the O(a^) full diagrams 
G^ N \ are irreducible. The 0(a s ) examples are Figs. 2(d) and 2(e). The factorization of reducible diagrams is the 
same as that of Figs. 2(a)-2(c): we simply insert the Fierz identity to separate G^ JV+1 - ) into the convolution of G^ N ' 
with the 0(a s ) collinear divergent diagrams, which contribute to the pion wave function cf)^>. The result is shown in 
Fig. 4(a). The explicit expression of has been derived in Eq. @. 

We then discuss factorization of irreducible G^ N+1 ^ . Let a be the outer most vertex on the upper quark line, and j3 
denote the attachments of the other end of the radiated gluon inside G^ N \ The pion possesses the leading-twist spin 
structure /Pi 75 from the Fierz identity. The fermion propagator adjacent to the vertex a is proportional to /Pi in the 
collinear region with the loop momentum I parallel to P\. Hence, the component 7+ in 7", which is located between 
the spin structure and the fermion propagator, gives the leading contribution. The vertex (3 must be dominated by 
the minus component. With the above reasoning, we propose the following replacement for the tensor g a @ appearing 
in the propagator of the radiated gluon, 

n (97) 
9af3 -> r ■ (27 

n • t 



The light-like vector n a in the minus direction, which was introduced in Eq. (|13|), selects the plus component of 7". In 
the collinear region I lies mainly in the plus direction, and lp selects the minus component of the vertex (3. Therefore, 
Eq. ( pjj ) extracts the leading-twist collinear divergences from the irreducible G^ N+1 \ 

The contraction of I 13 hints the application of the Ward identity in Eq. (^3|) to the case with two external on-shell 
quarks. Figure 4(b) contains a complete set of contractions of I 13 , which are represented by arrows, since the second 
and third diagrams have been added. The cuts on the quark lines denote the insertion of the Fierz identity. The 
Ward identity, namely, the equation described by Fig. 4(b), holds. The second diagram gives 

' [/ - (1 - x) A + (1 - x) /Pi] A 75 , 



(1-x) /Pi- / 
(1 - x)P 1 2 75 

(i-x) A- / 



A75 • (28) 



The first term in the second expression vanishes because of the on-shell condition Pf = 0. The second term corresponds 
to the 0(a^) full diagrams G^ N \ Similarly, the third diagram leads to 

where the first term vanishes and the second term corresponds to G^ N \ The factor n a /n ■ I from the collinear 
replacement in Eq. ( p7| ) is exactly the Feynman rule associated with the eikonal line in the direction of n. Equations 
( p8| ) and ( |29| ) imply that in the collinear region the irreducible G ( - N+1 ' ) are factorized as shown in Fig. 4(c). Obviously, 
the factorization of the irreducible diagrams with the gluon emitted from the outer most vertex on the lower quark 
line exists. We conclude that the irreducible G' Ar+1 ' can be factorized into the convolution of G^* 1 with described 
by Figs. 3(a)-3(d), in which the radiated gluon attaches the eikonal lines. 
We derive the collinear factorization of G^ N+1 \ 

G( N+ V w 4F> (8 GW , (30) 

where <p^ contains both the reducible gluons in Fig. 4(a) and the irreducible gluons in Figs. 3(a)-3(d). The remaining 
part F( JV+1 ), defined via 

G( W+1 ) = ^®GW + F( W ', (31) 

is infrared finite, i.e., free from the collinear divergence associated with the pion. The above procedures are also 
applicable to the 0(a^ +1 ) pion wave function cj>( N+1 \ which is defined by the perturbative expansion of Eq. (p^). 
We have 

f+^^gfl+F^ 1 ), (32) 



G 



with the infrared-finite function F^ N+1 \ 

Employing Eqs. (ES), (j31|), and (Bgh, we write 



AT 



+ _ ^(J+1)] g, + F (iV+l) 



J=0 



£ 0« ® i7( JV + 1 -J) + ff^+D , (33) 
j'=i 



with the 0(a^ +1 ) hard amplitude, 



AT 



H (N+i) = F (iv+i) _ ^ ^o+i) g, ff (Ar-j) _ (34) 

3=1 

Obviously, the function H( N+1 } does not contain any collinear divergence, since both F and F are infrared finite. 



Because of (j>(°> = 1, Eq. (33) becomes 



AT+l 

G (w+i) = J2 4>W ® i/( JV + 1 ^) . (35) 

The above expression indicates that all collinear divergences in the full diagrams of irj* — > 7 can be factored into the 
definition of the pion wave function in Eq. (^2|) order by order, and that the remaining hard amplitude is infrared 
finite. Allowing N to approach infinity, we prove factorization theorem for the process 777* — * 7 to all orders. 

At last, we prove by induction that soft divergences do not exist in the process 777* — > 7 to leading twist. The 
0(a s ) soft cancellation has been explained in Sec. II. Assume that the O(a^) full diagrams do not contain any 
soft divergence (though they contain collinear divergences). Consider the 0(a^ +1 ) full diagrams G^ N+1 K Similarly, 
we look for the gluon radiated from the outer most vertex on the upper quark line, and classify the diagrams into 
the reducible and irreducible types. For reducible diagrams, we express G( N+1 ) as the convolution of G^ with the 
0(a s ) infrared divergent diagrams as shown in Fig. 4(a). It is easy to confirm that the derivation in Eqs. (||)-(|8|) still 
applies to the soft region with the loop momentum I — > 0. The soft divergences in the 0{a s ) diagrams cancel in the 
same way as in Figs. 2(a)-2(c). The remaining 

G {N) 

are free from soft divergences as assumed above. Hence, the 

reducible G^ N+V> have no soft divergences. 

We then consider irreducible G^ N+1 \ The diagrams G^ , without soft divergences, are dominated either by hard 
or by collinear dynamics. In the hard region of G^- N \ internal particle propagators behave like 1/Q 2 . The attachment 
of a soft gluon, producing one more hard propagator, does not introduce soft divergences. The reason for this absence 
of soft divergences is the same as in Eq. (|To|). In the collinear region of G^ N \ where momenta parallel to Pi dominate, 
we employ the eikonal approximation for the small loop momentum I, 

The contraction of the numerator P\ a to the outer most vertex 7", which is mainly 7 + , leads to a vanishing contri- 
bution. Therefore, the irreducible G^ N+1 ^ do not contain soft divergences either. Extending N to infinity, we prove 
that the process irj* — > 7 is free from soft divergences. 



IV. FACTORIZATION OF tt 7 * -» n 



We investigate infrared divergences in the process irj* — > 71", and discuss only the factorization of the initial-state 
pion wave function. The discussion of the final-state wave function is the same. Assume that the incoming and 
outgoing pions carry the momenta Pi and P2, respectively, which are defined by Eq. (||), and that the momentum 
fraction Xi (^2) is associated with the lower quark line in the incoming (outgoing) pion. The above kinematic variables 
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have been indicated in the lowest-order diagrams in Fig. 5. Similarly, we consider the region with large momentum 
transfer Q 2 = —q 2 , q = P2 — Pi being the virtual photon momentum. 

Contracting the four-quark amplitudes in Fig. 5 with the leading-twist spin structures /Pi 75 J^/2N~ C and 75 JP<xl\/1N c 
from the initial and final states, respectively, we derive 

H 5a (xi,x 2 ) = -eg C F - _ _ ^ = ~^q^ . (37) 



„(o), x i 2n A7/x(/Pi - x 2 /P2)lu /Pi 75] 4eg 2 C F . 



The other two lowest-order diagrams, where the virtual photon attaches the lower quark lines, lead to the same 
expressions but with different electric charge e. Since contributions from soft partons (ari, X2 — * 0) are suppressed by 
the pion wave functions, the exchanged gluons in Fig. 5, off-shell by 0(Q 2 ), are regarded as being hard. 



A. 0(a s ) Factorization 

We first identify infrared divergences from 0(a s ) radiative corrections jl3| to Fig. 5(a). The diagrams in Fig. 6 
contain potential infrared divergences associated with the incoming pion. We do not consider self-energy corrections 
to the incrnal lines, since they, without infrared divergences, give only a next-to-leading-order hard amplitude. The 
treatment of Figs. 6(a)-6(c) is the same as that of Figs. 2(a)-2(c): soft divergences cancel and collinear divergences 
are factored into the pion wave function tp^ 1 ' by inserting the Fierz identity. 

Diagrams with the additional gluon attaching the internal lines, such as Figs. 6(d)-6(g), do not generate soft 
divergences, because the internal lines are off-shell. We concentrate only on collinear divergences. Following the 
detailed calculation in Appendix A, the collinear divergences in Figs. 6(d) and 6(e) are given by 



^ ^^Ma^l^-i^l,^)] , (39) 



he « ^^J[Hi°J(x 1 ,x 2 )-H^(^x 2 )] , (40) 

with £1 = x\ — l + 1 Pi . The dependences on l~ and It in have been neglected, since they are subleading in 
the collinear region. Collinear divergences from I parallel to Pi vanish in Figs. 6(f) and 6(g). We have separated 
the integrands I$d and I§ e into two terms, similar to those in Eqs. ( |l4|) and (|l7|). However, the color factors are 
different because of the different color flows in Figs. 2(d) and 6(d). We shall explain how to fix the color factors after 
calculating Figs. 6(h)-6(k). 

The collinear factorization of the irreducible corrections in Figs. 6(h)-6(k) relies on the cikonal approximation. The 
results are written as 

hh * ~2^n/^ H ^ (XuX2) ' (41) 
hi^^^/^H^^x,), (42) 

hj^^^^H^i^x,), (43) 

hk^-^^/^H^ixux,). (44) 

To discuss soft divergences, we simply allow l + to appraoch zero, i.e., £1 — ► X\. It is easy to find that soft divergences 
cancel between Figs. 6(h) and 6(i) and between Figs. 6(j) and 6(k). Hence, irreducible corrections are free from soft 
divergences. 

Combining Eq. (^|) and Eq. (|4|) with the first and second terms in Eq. (^9|) , respectively, and combining Eq. ( |4^ ) 
and Eq. ( j44| ) with the second and first terms in Eq. (p0|), respectively, we obtain the correct color factors: 

I 6 (d)- m « <t$[H&\xi,X2) - Hgfaxu)] 

+4>i X hHf]{x 1 ,x 2 )-H^{i 1 ,x 2 )}, (45) 
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which are also described by Figs. 3(a)-3(d). It is observed from the above expression that the collinear divergences in 
the summation of Figs. 6(d)-6(k) are identical to those in the summation of Figs. 2(d) and 2(e), consistent with the 
universality of hadron wave functions. 

We then investigate infrared divergences from radiative corrections to Fig. 5(b). The diagrams in Fig. 7 contain 
potential infrared divergences associated with the incoming pion. The factorization of Figs. 7(a)-7(c) is the same 
as that of Figs. 6(a)-6(c): They are absorbed into the leading-twist pion wave function (fP^ straightforwardly. The 
contributions from Figs. 7(d) and 7(e) are the same as those from Figs. 6(d) and 6(e). 

For Figs. 7(f), we consider only the collinear divergences from I parallel to Pi, 

T ~ ig2 n-[(l-a;i)Pi+/] (0) 

l7f ~ N- c [(l-x 1 )P 1 +l]Wn.(P 1 +l) H5b {XUX2) ■ (46) 

The result differs form that of Fig. 6(f), which vanishes. Note that the collinear divergence in Eq. ( fili| ) is not in the 
correct eikonal form. To show that it does not cause any trouble, we consider the collinear factorization of Fig. 7(h), 

t r^W 2 n-[{\- x l )P l +l] re-fi (Q) 

l7h " ^ [(1-^)^+^.(^+0 n .l H »> ^'^ ■ (4?) 
Combining Eqs. (^6|) and (E7l), we have the simpler expression 

hf+7h * -2^n/^ H ^ [XuX2) ' (48) 

whose collinear divergence is exactly the same as that of Fig. 6(h) in Eq. pi]). At last, combining the collinear 
divergences in Eq. ( |4g| ) and those from Fig. 7(d), we obtain the result described by Fig. 3(a). 

A similar analysis shows that the collinear divergences in the sum of Figs. 7(g) and 7(j) are the same as in Fig. 6(j). 
Combined with the divergences from Fig. 7(e), we obtain the result described by Fig. 3(c). The results for Figs. 7(i) and 
7(k) are the same as for Figs. 6(i) and 6(k), respectively, and their collinear factorizations, combined with Figs. 7(d) 
and 7(e), are described by Figs. 3(b) and 3(d), respectively. We conclude that the 0(a s ) collinear divergences in 
Figs. 6 and 7 are identical to those in Fig. 2, and can be formulated into the leading-twist pion wave function cfft' 
defined in Eq. ©. 



B. All-order Factorization 



The analysis of the 0(a s ) collinear divergences in ttj* — > it indicates that factorization of a QCD process requires 
summations of many diagrams. For example, the summation of Figs. 7(f) and 7(h) removes the dependence of the 
hard amplitude on the loop momentum I, The summations of Figs. 6(d) and 6(h) and of Figs. 7(d), 7(f) and 7(h) 
give the correct color factors. To prove factorization theorem to all orders, we deal with these summations by means 
of the Ward identity. It is trivial to generalize the proof for the process ir^* — > 7 to tt^* — ► w. In the present case the 
amplitude G M in Eq. ( p3| ) contains four external on-shell quarks. We simply interpret the function H as the part that 
does not involve the collinear divergences associated with the incoming pion. Then repeat the procedures and further 
factorize H into the the convolution of a hard amplitude with the wave function associated with the outgoing pion. 

The proof of the absence of soft divergences in the process irj* — > it is subtler. It has been shown that the 0(a s ) 
corrections do not produce soft divergences. Assume that the soft cancellation exists up to 0(a^). For reducible 
0(a^ +1 ) full diagrams G^- N+1 \ we insert the Fierz identity to factor out the 0(a s ) infrared divergent diagrams, 
which involve the gluons emitted from the outer most vertices on the upper and lower quark lines. This insertion 
works for both the initial- and final-state pions. The soft divergences in the 0{a s ) diagrams cancel in the same way 
as in Figs. 2(a)-2(c). The remaining O(a^) full diagrams G^ have no soft divergences as assumed above. Hence, 
the reducible G ( - JV+1 - ) do not contain soft divergences. 

For irreducible G l -- N+1 \ we consider the gluon emitted from the outer most vertices on the upper and lower quark 
lines in the incoming pion. Since 

G (N) 

are free from soft divergences, they are dominated either by hard or by collinear 
dynamics. In the hard region and in the collinear region with momenta parallel to P±, the soft divergences vanish 
for the same reason as for the process ttj* — > 7. In the collinear region with momenta parallel to P2, we adopt the 
eikonal approximation, 

#^7. p2~—,p2- (49) 
(Pi + i) n ■ I 
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It is easy to observe that the above approximation can be achieved by the replacement in Eq. (|27|). 

We propose to employ Eq. (|27| ) to extract the soft divergences in the irreducible G^ N+1 \ In the hard region and in 
the collinear region with momenta parallel to Pi , the replacement in Eq. ( p7| ) simply modifies subleading contributions, 
that are free from soft divergences. In the collinear region with momenta parallel to P2, the replacement extracts 
the correct soft divergences in the irreducible G^ N+1 K Hence, Eq. (|27j) always works for a leading-twist analysis 
of soft divergences. The derivation from Eq. © to Eq. © holds, and the irreducible G^ 1 ' are factorized into 
the convolution of G^ with Figs. 3(a)-3(d). Since G^ have no soft divergences, we allow the loop momentum I 
flowing inside to approach zero. It is then obvious that the soft divergences cancel between Figs. 3(a) and 3(b) 

and bwteen Figs. 3(c) and 3(d). That is, the irreducible G^ N+1 \ like the reducible G^ N+1 \ are also free from soft 
divergences. The same argument applies to the gluons emitted from the outgoing pion side. At last, extending N to 
infinity, we prove the absence of soft divergences in the process 777* — ► tt. 

We conclude from the discussions in Sees. Ill and IV that the leading-twist infrared structures in the processes 
777* — ► 7r and irj* — > 7 are identical, and that the pion wave function, defined by Eq. (123), is universal. 



V. FACTORIZATION OF B MESON DECAYS 



In this section we present the all-order proof of factorization theorem for the exclusive B meson decays B — > 
and B — > ttW in the heavy quark limit. It will be shown that the proof for the pion transition form factor in Sec. Ill 
can be generalized to B meson decays, if terms suppressed by powers of 1/Mg are neglected. The momentum Pi of 
the B meson and the momentum Pi of the outgoing on-shell photon (pion) are parametrized as 

Mf> , * „ Mr , , 

P 1 = - 7 |(1,1,0 T ), ft=-^(0,Tj,O r ), (50) 

where r\ denotes the energy fraction carried by the photon (pion). Assume that the light spectator quark in the B 
meson carries the momentum k, and that e denotes the polarization vector of the photon. We consider the kinematic 
region with small q 2 , q = Pi — P2 being the lepton pair momentum, i.e., with large 77, where PQCD is applicable. In 
the heavy quark limit the mass difference between the B meson and the b quark, A = Mb — rrib, is a small scale, which 
will appear as higher- twist terms proportional to A/ Mb- The four components of the spectator quark momentum k 
are of the same order as A. 

The lowest-order diagrams for the B — > decay and for the B — ► irlv decay are displayed in Fig. 1 and Fig. 5, 
respectively, but with the upper quark line replaced by a b quark line and with the vertex x representing a weak 
decay vertex. It is easy to observe that Figs. 1(a) and 1(b) scale like 

H '°' " 5 - k) 2 - Ml ~ Ml • (52) 

indicating that Fig. 1(b) is power-suppressed. In the following analysis we shall concentrate only on Fig. 1(a). There 
are two leading-twist spin structures for the B meson, 

( /Pi + M B ) A75 {JPx + M B ) Als 



2JN r 2JN r 



(53) 



with the dimensionless vecor n = (1,0, 0t), each of which is associated with a B meson wave function. For the 
radiative decay B — > jlP, only the first structure contributes at the leading- twist. 
Contracting the amplitudes in Fig. 1(a) with (/Pi + M b) A75 / '(2\fW c ) , we derive 

r (0), s _ VN~ctr[jt /P 2 7m(1-75) A] 



where x — k + / P^~ is the momentum fraction. Since P2 has been chosen in the minus direction, only the plus 
component k + of k is relevant. We have dropped the higher- twist terms proportional to k in the numerators and A 
by assuming mj « Mb- It is then possible to integrate out the k~ and fcx dependences in the B meson wave function 
■0+(fc), leading to the light-cone B meson wave function, 

<p+(x) = J dk~d 2 k T ^+(k) . (55) 

It will be shown that a light-cone B meson wave function can always be defined, if an appropriate frame is chosen, 
even though the four components of k are of the same order. 
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A. 0(a s ) Factorization of B — > ■ylu 

We start with the one-loop diagrams in Fig. 2 for the B — > 7/1/ decay, discussing the factorization of their infrared 
divergences. The factorization of Figs. 2(a)-2(c) requires the insertion of the Fierz identity. Take Fig. 2(b) as an 
example (the analysis of Figs. 2(a) and 2(c) is trivial). The loop integrand is written as 

T ~ VwJ / 757° P\- fc+ /+M B v (J> + I)frft\ l 

hb ~ 19 CFtr y»-[k=w— ( Pl - k + if - Mf — 2 — )p 

X€ ~2 (P 2 - k + 2 ' ( j 

where we have kept only the leading-twist structure 757° in the Fierz identity, and introduced the dimensionlcss 
vector v = Pi /Mb- The other structures, such as 75 and 75(7"^, do not contribute, because the corresponding hard 
amplitude vanishes: 

^VKtrlMA- /)7m(! - 7 5 ) M B75 ] _ Q 



2 (P 2 - k + I) 2 

In the heavy quark limit we have the eikonal approximation, 

/Pi- /&+ / + Mb „, , , ~ v 1 1 i jl r\ t*st\ 
(fi-fc + O'-Ml 7 ( ^ + 7) * MT^fc) W + 7 ) > ( 58 ) 

which is in fact equivalent to the heavy quark expansion. Obviously, only the component 7+ = 7~ of 7 Q contributes 
to the hard amplitude at leading twist, requiring the spin structure 757° to be 757 + in the first trace for the B meson 
wave function. Equation ( |56| ) becomes 

h^i^UH^iO, (59) 
with £ = x — l + 1 'Pi and the 0(a s ) B meson wave function, 

(0+ U = *S C F trl Ivjj—^-z j h^T^k) ■ (60) 



Figure 2(b) has been expressed as the convolution of the hard amplitude with the 0(a s ) infrared divergent 

diagrams, which contribute to <f>v , in the momentum fraction £ 
Equation ( |60| ) is simplified into 



Performing the contour integration over, say, l~ , we observe that the integral is singular only when the component l + 
is of O(A). This observation implies that the infrared divergence associated with the B meson is of the soft type. This 
dynamics differs from the collinear type of divergences associated with the pion in the process 717* — > 7. It is easy to 
understand that the soft divergences in Figs. 2(a)-2(c) do not cancel in B meson decays 0: The light spectator quark, 
carrying a small amount of momenta, forms a color cloud around the b quark. This cloud is also huge in space-time, 
such that soft gluons resolve the color structure of the B meson. 

Diagrams with the additional gluon attaching the internal quark in Figs. 2(d) and 2(e) also contain soft divergences, 
since the internal quark is off-shell only by O(AMb)- Note that the internal quark in the process ir-f* — > 7 is off-shell 
by 0(Q 2 ). There is no collinear divergence from I parallel to P 2 , because the internal quark propagator remains the 
scaling law 1/AMb, instead of 1/A 2 . The integrand associated with Fig. 2(d) is written as, 

h -'= -'"■><>■— tr \ ^jp—w^ { P 2 -k + i)^ {l - l5) 
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Inserting the Fierz identity, we have, 



l2d 



-ig 2 C F 



2P 2 u 



. I 757" fa~ h± [±Mb V (fi + I) fa I 1 



(P 2 -£; + 2 2 {Pi-k + l) 2 ~ M 



i 2 



xe- 



fW c tr[jt /P 2 7m( 1 ~ 75)M B7a 7 5 ] 
2 (P 2 - A;) 2 



(63) 



Obviously, the first trace is proportional to n a in the heavy quark limit. 
We employ the relation similar to Eq. (p~3|), 



(P 2 -k + l) 2 n-l 



1 



(P 2 - fc) 2 



(p 2 - * + 



(64) 



where the terms (k — l) 2 ~ (3(A 2 ) and /c 2 ~ <3(A 2 ) have been neglected. The eikonal line in the direction of n appears. 
Equation (63) reduces to 



with the 0(a s ) B meson wave function, 

757 + + 1 1 



■ + >hd = -ig'C F -^- l tr< 



Pv-(l-k) 



(65) 



(66) 



The above expression implies that the infrared divergences in irreducible diagrams can also be collected by the eikonal 
line along the light cone. This is attributed to the choice of the frame, in which the photon moves in the minus 
direction. The first and second terms in Eq. (|6^) are described by Figs. 8(a) and 8(b), respectively, where the eikonal 
lines in the directions of v and of n have been indicated. In Fig. 8(b) the loop momentum I flows into the internal 

quark line, such that the second term appears as a convolution of with (4>^) 2 d in the momentum fraction £. 



Figure 2(e) gives the loop integrand 



/Nr. 



he = ieg C F ——tr< 7 



v Jk- J, fa- fc+ / 



2 

fa- h 

(P2 - k) 2 



4 



(fc-0 2 ^fa-k + iy 



7^(1™ 75) (/Pi +M B )Al5 



(67) 



Inserting the Fierz identity, we have, 



he ~ ig Cf 



2P 



2v 



{p 2 -k + iy 



:tr< 7 



v tl I 757 a (/' + ^)A75 \ 1 
(k-l) 2 2 2 



/ 2 



xe- 



^K tr\J fa ltl {\ - 75)M B7a 7 5 ] 
2 (P 2 - k) 2 



(68) 



Similarly, the component 7+ = 7 of j a in the second trace is selected. Employing Eq. (p4|), the above expression is 
rewritten as 



h 



with the 0(a s ) B meson wave function 



<^)2e[H^\0~H^( X )], 



he = ig 2 C F ^tr< 7 



n ■ I 



„ h- / 757 + (^ + PjA75 1 
(k-l) 2 2 2 [ I 2 



(69) 



(70) 



We have split Fig. 2(e) into two terms, whose infrared divergent factors are described by Figs. 8(c) and 8(d), respec- 
tively 
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Comparing Eqs. (^0|), (^6|) and j7C|), the Feynman rules for the perturbative evaluation of the B meson wave function 
are clear: cj)g can be written as a nonlocal hadronic matrix element with the structure 7s7 + /2 sandwiched: 

$\x) = J %^e" p iy- myn^i-ig) J" dzn-A{zn)b v mB{Pi)) , (71) 

where the rescaled b quark field, 

b v (y) = exp(iM B v ■ y)^^b{y) , (72) 

has been introduced. The Feynman rules associated with b v are those for an eikonal line in the direction of v given 
in Eq. (|58|). The above definition reproduces the contributions from Figs. 2(a)-2(c) and from Figs. 8(a)-8(d), if it is 
evaluated perturbatively. 

The factorization formula for the B — > decay is written, up to 0(a s ), as 

(0f + ® (H[°J + H$) + 0{al) , (73) 

with <fi\/ = 1 and ® representing the convolution in the momentum fraction. The 0(a s ) hard amplitude H^J is 
defined by, 

where the subscript i runs from 2a to 2e. Obviously, if}* is infrared finite, since all the (a s ) soft divergences have 
been absorbed into the B meson wave function (frP . 



B. All-order Factorization of B — > *ylv 



We prove leading-twist factorization theorem for the B — ► jlv decay to all orders, and construct a gauge-invariant 
light-cone B meson wave function defined by 



"- <•<■> 7^T7 / ^e^ + ^(0| g -(^)^Pexp 



/ cfen • A(zn) 



fc„(0)|B(Pi)) , 

(75) 



as shown in Fig. 8(e). The eikonal lines in the directions of v and of n have been indicated. By expanding the light 
quark field q(y~) and the path-ordered exponential into powers of y~ , the nonlocal matrix element can be expressed 
as a series of covariant derivatives (D + ) n q(0), implying that Eq. ([75]) is gauge invariant. 

We present the proof by induction. The factorization of the 0(a s ) infrared divergences associated with the B meson 
has been worked out. Assume that the factorization of the infrared divergences holds up to O(a^), that is, we have 
Eqs. (E^)-(E^) for the B — > decay. Consider a complete set of 0(a^ +1 ) full diagrams G^ N+1 K We look for the 
gluon, one of whose ends attaches the outer most vertex on the b quark line. We classify G' Ar+1 - ) into the reducible 
and irreducible types as in Sec. III. The factorization of the reducible G^ N+1 " > is the same as that of Figs. 2(a)-2(c): 
Following Eqs. (|5^)-(|60|), we insert the Fierz identity to separate the reducible G^ N+1 ^ into the convolution of 

with the 0(a s ) infrared divergent diagrams, which contribute to the B meson wave function (f>^ ■ The result is similar 
to that shown in Fig. 4(a). 

For the factorization of the irreducible diagrams, we rely on the Ward identity in Eq. (^3j). In this case the amplitude 
G M contains two on-shell external quarks, one of which is the heavy b quark. As hinted by the 0(a s ) analysis, the 
soft divergence, associated with the gluon radiated by the outer most vertex on the b quark line, can be collected by 
the eikonal approximation, i.e., by the replacement in Eq. (|2~7j). The reason is as follows. In the heavy quark limit 
the full diagrams of the B — > ^Iv decay are dominated by the momentum flow along the photon momentum P% in 
the minus direction. Strickly speaking, they are dominated by Pi — k. Hence, the vertex /3 inside G^ the radiated 
gluon attaches is mainly minus, and the vertex a on the b quark line is mainly plus. Apparently, the tensor n a lp/n ■ I 
extracts the correct leading contribution. Therefore, we have the equation described by Fig. 4(b). 

The second and third diagrams in Fig. 4(b) give 
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Pi- jk- /- M B 

= -~B , 1 , M M/> + I) fil* -(/O + I) /175 , (76) 

/Pi- jk- J, -M B 

l p (/b + I)jhrft>f-L- = (jO + I)fi l5 yft-r^-^ + ^/iTfe, (77) 

respectively. After assigning the quark propagators on the right-hand sides of the above expressions into the corre- 
sponding loop integrals, the first terms are proportional to k ~ 0(A). They are suppressed by 0(A/Mb) compared to 
the second terms, which correspond to G^ N \ Neglecting the higher-twist terms, the irreducible G^ N+1 " > are factorized 
into the convolution of GW with the 0(a s ) B meson wave function, in which the radiated gluon attaches the eikonal 
lines as in Figs. 8(a) and 8(b). The factorization of the irreducible G^ N+1 \ with the gluon emitted from the outer 
most vertex on the light spectator quark line, is similar. The resultant 0(a s ) infrared divergent diagrams for the B 
meson wave function are those in Figs. 8(c) and 8(d). 

Combining the factorizations of the reducible and irreducible G^ N+1 \ we arrive at Eq. ([3lf). The factorization in 
Eq. (32) for the B meson wave function 0^ v+1 - ) also exists. Following the steps in Eqs. (|33|)-(|35|), we complete the 



all-order proof of leading-twist factorization theorem for the B — > jlD decay. The definition of the hard amplitude 
at each order is the same as in Eq. (|34|). Since the full diagrams are dominated by momenta along P2, only the plus 
component k + of k is relevant in the hard amplitude. This is the reason we can integrate the B meson wave function 
over k~ and fey, obtaining the light-cone B meson wave function in Eq. (f75|). 



C. Factorization of B — » nil/ 



The factorization of the pion wave function in Sees. II, III, IV, and the factorization of the B meson wave function 
in the B — > decay can be applied to the B — > ttW decay straightforwardly. Here we simply explain some points 
of theproof, and neglect the details. The lowest-order diagrams are shown in Fig. 5. Both the spin structures in 
Eq. ( p3| ) contribute to Fig. 5(b), but only the first one contributes to Fig. 5(a). Since Pi has been chosen in the 
minus direction, only the plus component k + of k is relevant. It is then possible to define a light-cone B meson wave 
function for the B — ► nlv decay. We first identify infrared divergences from 0(a s ) radiative corrections to Fig. 5(a), 
which are displayed in Fig. 6, since their analysis is similar to that of the B — > 7ZP decay. Note that Fig. 6 is complete 
only for the construction of the B meson wave function. We do not consider self-energy corrections to the internal 
lines, which give only next-to-leading-order hard amplitudes. 

The construction of the pion wave function is basically the same as that for the process ir-f* — > ir. That is, the 
infrared divergences associated with the outgoing pion is of the collinear type. Soft divergences cancel by pairs, 
because soft gluons do not interact with the color singlet pion. For example, soft divergences cancel between Fig. 6(h) 
and 6(i) and between 6(j) and 6(k). The collinear gluons are still collected by the eikonal lines along the light cone, 
even though they may attach the heavy b quark. The reason is that when a loop momentum I is parallel to the pion 
momentum P2 in the minus direction, only the component 7+ of the vertex on the b quark line and P x + of the B 
meson momentum are selected. Consequently, we derive the definition for the pion wave function, which is identical 
to Eq. (|2^). This conclusion is consistent with the universality of hadron wave functions. 

We then concentrate on the B meson wave function. The treatment of the reducible diagrams in Figs. 6(a)-6(c) is 
exactly the same as of Figs. 2(a)-2(c) for the B — > jli> decay. We insert the Fierz identity to separate these diagrams 
into the convolution of Fig. 5(a) with the 0(a s ) soft divergent diagrams, which contribute to the B meson wave 

function <jf^ . The investigation of the irreducible diagrams in Figs. 6(d)-6(e) hints that the soft divergences can be 
collected by the eikonal lines along the direction of n, a conclusion similar to that for Figs. 2(d) and 2(e). Hence, the 
replacement in Eq. (|27|) for the loop gluon extracts the soft divergences of the irreducible diagrams in Figs. 6(d)-6(k). 
The Ward identity applies, and the sum of these irreducible diagrams is factorized into the convolution of Fig. 5(a) 

with part of 4>^ described by Figs. 8(a)-8(d). Following the procedures in Sec. Ill, we prove leading-twist factorization 
theorem for Fig. 5(a) in the B — > ttW decay. The B meson wave functions 4> + constructed from the B — > jlv and 
B — > ttW decays are identical. 

The discussion of the 0{a s ) radiative corrections to Fig. 5(b) shown in Fig. 7 is similar, though more complicated. 
The soft approximation, the Fierz insertion, and the Ward identity apply. Because of the existence of the heavy quark 
propagator, the spin structures 757°, 75 and 750""^ in the Fierz identity contribute. In this case two leading- twist B 
meson wave functions can be constructed Besides <fi + defined in Eq. ((7^), we derive the additional light-cone B 
meson wave function associated with the structure (f \ + Mb) jh"f^/(2y/Nc), 
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M*) 



^e^ + ^(0| g -(^)^-Pexp 

Z7T Z 



-z<7 ; dzn ■ A(zn) 
'o 



MO)|B(Pi)) . 



(78) 



In summary, the nonlocal hadronic matrix element for the B meson is expressed as 



^e" p 1 + ^(0|g 7 ( y -)6M0)|B( J Pi)) = 



75 



(79) 



7/3 



This result is exactly the same as that obtained in || at leading twist. 

A remark is in order. The hard amplitude from Fig. 5(b), proportional to X2 being the momentum fraction 

associated with the pion, develops an infrared singularity, if the pion wave function vanishes like X2 as x-i —> [ fTI] . 
Even though we have proved leading- twist factorization theorem for the B — > nlv decay, its practical application is 
questionable. It has been shown that the inclusion of parton transverse momenta kx smears the singularities from the 
end-point region of momentum fractions @] . When kx is included, the dependence on the transverse loop momentum 
It, being of the same order as kr, is not negligible in the hard amplitude. This complexity can be resolved by Fourier 
transforming the factorization formula from the kx space into the b space, where b is the impact parameter conjugate 
to kr (l^|l(|. The It dependence then appears in the factor exp(— Ht ■ b), and decouples from the hard ampltude. 
This factor can be assigned into the corresponding loop integral, which contributes to the definition of the meson 
wave functions. We claim that the factorization of the B — > irW decay must be performed in the impact parameter 
space fTl. 



VI. CONCLUSION 



In this paper we have investigated the infrared divergences in the processes ttj* — > 7 and irj* — ► n. We summarize 
their comparision below. There are no soft divergences, since they are either absent or cancel among sets of diagrams. 
In the collinear region with I parallel to Pi, Figs. 2(a)-2(c) are identical to Figs. 6(a)-6(c) [Figs. 7(a)-7(c)]. Figures 2(d) 
and 2(e) are identical to the combination of Figs. 6(d)-6(k) [Figs. 7(d)-7(k)]. That is, the collinear structures are the 
same at the leading twist, consistent with the concept of universality of hadron wave functions in PQCD factorization 
theorem. However, due to the potential significant sublcading contributions at low energies, the extraction of the 
leading-twist pion wave function from the ttj* — > j(tt) data suffers ambiguity. For details, refer to Appendix B. 

We have presented an all-order proof of leading-twist factorization theorem for the processes ttj* — > 7 and irj* — ► n, 
and for the decays B — > jlD and B — > nlv by means of the Ward identity. The small scales, such as the light spectator 
quark momentum k and the B meson and b quark mass difference A, are neglected in the heavy quark limit, such that 
the Ward identity holds. Our proof is simple compared to that in 0], and explicitly gauge invariant, compared to that 
in [Q. We have constructed the gauge-invariant pion and B meson wave functions, and confirmed their universality. 
The path-ordered integral appears as a consequence of the Ward idenity. Note the difference between the definitions 
for the pion and B meson wave functions. The former, collecting the collinear divergences, depends on the structure 
757 + . The latter, collecting the soft divergences, depends on the structure 75 /v. We emphasize that it is possible 
to define a light-cone B meson wave function, if an appropriate frame is chosen, in which the photon (pion) moves 
in the minus or plus direction. This is also the reason we can extract the infrared divergences in n^f* — > 7(71") and 
B — > -f(ir)lD using the same replacement in Eq. (p7j). 

The leading-twist factorization of the B — > nlv decay can be proved straightforwardly, following the procedures 
presented in Sees. II-V. However, for a practical application, the parton transverse momenta kT must be included 
in order to smear the end-point singularities in the hard amplitude. We shall derive factorization theorem including 
these additional degrees of freedom elsewhere Our proof will be generalized to nonleptonic B meson decays, such 
as B —> 7T7T. This factorization is more complicated, since nonleptonic decays involve three characteristic scales: the 
W boson mass Myy, Mb, and small scales of 0(A), such as /ct- 
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discussions. This work was supported in part by the National Science Council of R.O.C. under the Grant No. NSC- 
89-2112-M-006-033, by Grant-in Aid for Special Project Research (Physics of CP Violation) from the Ministry of 
Education, Science and Culture and by Hiroshima University, Japan. 
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APPENDIX A: 0(a s ) COLLINEAR CORRECTIONS 



In this Appendix we supply the details of the derivation of the collinear divergences in the process ttj* — * n. The 
loop integrand from Fig. 6(d) is written as 

T ~ %egA + S A p fl ft-Hl A± I (l-^l) A + / 

2A C ^ (P 2 - xiPi+/j p [(1 - xijPi + /J 2 

1 tr{T a T b T c )T ab S x 

X7 a /Pi75>7^ — p 5 n2 , p p-Ta, (80) 

i 2 (xiPi - a; 2 P 2 - J) ( x i "1 ~ 



with iV c being the number of colors and the triple-gluon vertex 

^abc 
a8\ 



K% = -f abc [g a0 (2l - X1P1 + x 2 P 2 ) x + g f3 x(2x 1 P 1 - 2x 2 P 2 - l) c 



+g\ a {x 2 P 2 -x l P l -l) ], (81) 

jabc j-jgjjjg a antisymmetric tensor. In the collinear region with I parallel to Pi, only the term proportional to gs\ is 
important, since the terms proportional to g\ a and g a 8) giving 

[(1 - u) A + Al a /Pi757 A .9A Q = 2Pi • [(1 - x l )P 1 + l] l5 ~ 0(A 2 ) , 

7" A7s(2 /-xi jPx+x 2 Abs A7 /? <?a/3 = 2Pi • (21 - x x P x ) /P 2 ~ 0(\ 2 ) /P 2 , (82) 

do not produce collinear divergences. Considering the second term in Eq. (|8l|), Eq. ( p0|) exhibits the collinear 
divergence 



-3 ?g 2 p 2 ■ [(l - x)Px + *] g(0) 

[(l-o;)Pi+q 2 Z 2 [P 2 -(Z-a;Pi)] J 



r j 2 [V -L - 1 -r ij W , , 

7 6d ~ 77 1 ;15>;2rn 71 Z77vt H 5a > I 8 '*) 



where we have employed the identities, 



tr(T a T b T c ) = ^(d abc + if abc ) , d abc / afcc = , / afcc / Qbc = 24 , (84) 

<P hc being a symmetric tensor, and the approximation for the denominator 

(a-iPi - x 2 P 2 - If = 2x 2 P 2 ■ (I - xiPi) . (85) 

It is found that the infrared divergent piece of the radiative correction has been completely factored out. The remaining 
it, denoted by the lowest-order hard 
Figure 6(e) gives the loop integrand 

T ie 9 4 , j a gl A~ / A d 3 A-_gl A 

/6« = 7T7T ir i 7 7 — 75 7717 75 A7 P 77i 7JT2 

2iV c ^ (xiPi-iy (P 2 -xiPi) 2 

] tr{T a T b T c )T ab S x 

X7mA75>727 — 75 75 7727 — 75 7772", (86) 

\l z (x 1 Pi-x 2 P 2 -iy(xiP 1 -x 2 P 2 y 



part, denoted by the lowest-order hard amplitude , does not depend on the loop momentum I at all 



with the triple-gluon vertex 

*/3A 



K% = -f abc [9aa{l - x x P x + x 2 P 2 ) x + g sx {2x x P x - 2x 2 P 2 + l) a 



+g\ a (x 2 P 2 - x x P x - 2l) p ) . (87) 
Similarly, in the region with I parallel to Pi we have 



T - 3t 9 o-(Q) 

6e ~ (nPj-i) 2 ! 2 5a ■ 



(88) 



The integrand associated with Fig. 6(f) is written as, 
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T 2£^1 + J t> c ^2~X 1 jf\ P 2 - X x P 1 + / 



X [(1 - xl)Pxt If Y A75 f Z 2 (xiPi - * 2 P 2 ) 2 ' (89) 



which is simplified into 



6/ ~ 9 F [(i - x l )p l + inp 2 ~ Xl p, + lyp^p, - x 2 p 2 r ■ 1 ' 

The trace in the above expression is proportional to the vanishing factor Pi - 1, which suppresses the divergence from 
the denominator [(1 — xi)P\ + I} 2 . Figure 6(g) does not contain collinear divergences for the same reason. 
The loop integrand associated with Fig. 6(h) is writen as 

Ieh = TXT — * r S 7a 75 Pilvjf. , p : , 12 7 jt, p , n2 7^ 

AN C y [(1 - a3 2 )P2 + if (P 2 - x%Pi + iy 



[(1 - a;i)Pi + Z] 2 ' ' '° - x 2 P 2 ) 2 ' 



which is simplified into 



2eg 4 C F P! ■ [(l-z 2 )P 2 +^ J A-.Ti^+Z 

6,1 _ w c [(i - x 2 )p 2 + i} 2 \ n (p 2 - + o 27 ' 1 

(i-^x) i f92) 

[(i-x 1 )p 1 + /] 2 J; 2 (a ;i p 1 -^ 2 p 2 ) 2 ■ iy j 

In the region with / parallel to P\, we have the approximation, 

j _»ff 2 P 2 -[(l-a:i)Pi+/] rr(0)^».9 2 rc-[(l-a:i)Pi+i] (0) 

6a ~ iV c [(1 - a^Pi + Z] 2 Z 2 P 2 • Z 4a 7V C [(l-xi)Pi+Z] 2 Z 2 n-Z 5a ' 1 J 

where we have adopted the approximation for the quark propagator [(1 — a; 2 )P 2 + I] « 2(1 — a; 2 )P 2 • /, and dropped 
all the terms proportional to P\ ■ I. 

The loop integrand corresponding to Fig. 6(i) is written as 

t e 9 A ° F + J X2 ^ 2+ I <* fa-xi JPi+ j 

/K = -W tr \ 7 " (x 2 P 2 + 2 >75 ^ 7 (P 2 -x lPl+ l)^ 



[(1 - .x^Pi + Z] 2 ' ' '° j PCuPi - x 2 P 2 - i) 2 ' 

whose sign is opposite to that of Eq. ([n|) due to the anti-quark propagator. The above expression is simplified into 

_ 2eg i C F \ x 2 P 2 + / P 2 - xj. f\+ / 
N c \{x 2 P 2 + l) 2 {P 2 -x 1 P 1+ l) 21 ^ 

It is easy to derive the collinear approximation, 



-ig 2 «■[(!- a:i)Pi + 1] {xiPi - x 2 P 2 f w (o) 
iV c [(1 - xi)Pi + l] 2 l 2 n ■ I (xiPi - x 2 P 2 - If 



l6l ~ AT \C\ _ ~^P.. J. 712/2^ . / („. R. _ „„p„ _ n2 H Sa ■ ( 96 ) 



The similar procedures apply to Figs. 6(j) and 6(k). 
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The integrand associated with Fig. 7(f) is written as 



which is simplified into 



l7f = -^ tr y a15 ^"(ft-w^^-^+o 37 

[(i - xi)Pi + i\ z i v{x\Pi — X2P2) 



2eg i C F P 1 -{l-x 2 P 2 ) J A - ^2 A 

7/ ~ iv c (Pi - x 2 p 2 + 2 \ (Pi x 2 p 2 y 

(l-ari) /Pi+ /l 1 



In the collinear region we have the approximation, 



The integrand corresponding to Fig. 7(h) is written as 



which is simplified into 



[(1 - x^Pt + if ' " ( p{ Xl i\ - x 2 p 2 y ' 



2eg i C F P 2 -[{l-x 1 )P 1 +l] j (1 - x 2 ) P 2 + / 
7 " A^c [(l-^Pi+Z] 2 r \[{l- X2 )p 2+ l]2^ 



We have the collinear approximation, 



(98) 



P 2 -[(l-xi)Pi+f] (o, (99) 



hh = ~^vT^\ 7q75 ^ [(l-x 2 )P 2 + l]^(P^x 2 P 2 + l)^ 



/Pi - X2 /p 2 + / pi 1 nrm 
x (Pi - x 2 p 2 + 2 J * 2 (*iPi - ^p 2 ) 2 ' 1 j 



j „ W 2 P2 ■ [(1 - zi)P + j] Pi ■ P2 (o) , . 

iv c [(i-xi)Pi + ;] 2 / 2 p 2 -(Pi + o p 2 -^ 56 ' 



APPENDIX B: DETERMINATION OF THE PION WAVE FUNCTION 



In this Appendix we comment on the determination of the leading-twist pion wave function, which can be 
parametrized as 



<fi(x) 



3/. 
JW, 



x(l — x) 



l + -c(5(l-2z) 2 -l) 



(103) 



with the shape parameter c. It has been proposed to extract the leading-twist pion wave function from experimental 
data of the pion transition form factor p8[. The asymptotic model was obtained with the shape parameter c ~ 0, 



x{\ — x) 



(104) 
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The pion wave function can also be determined from other processes involving pions, such as the pion form factor 
and the B meson decays B —* irlv and B —* Dir. It has been known that a large value of c is preferred for explaining 
the data of the pion form factor Q . 

Another quantity that has been considered is the ratio of the branching ratios of the B — > Dir dceays fp]^(| , 

„ B(B+ -> D°ir+) , , 

R ^ BiB^D-J) - < 105 > 

The charged B meson decay contains both factorizable and nonfactorizable amplitudes: the B — > D form factor 
associated with the external-VF emission, the B — > ir form factor associated with the internal-VF emission, and the 
nonfactorizable amplitude associated with the internal-IF emission. The neutral B meson decay contains only the 
factorizable external-VF emission amplitude. The Wilson coefficient for the factorizable internal- W emission amplitude 
is small at scales around the b quark mass. Hence, the difference between the branching ratios B(B + — > D°tt + ) and 
B(B° — > D~ir + ) is attributed to the nonfactorizable internal-VF emission amplitude. To explain the data of R ~ 1.6, 
a larger c ~ 0.5 has been obtained This value of c is located between those for the asymptotic model and for the 
Chernyak-Zhinitsky (CZ) model corresponding to c = 2/3 J22|, 

^{ X ) = ^k x (i- x )(i- 2x ) 2 . (106) 

Note that the coefficient c = 0.44 (at the factorization scale about 1 GeV) derived from QCD sum rules |2j| is close 
to that extracted from the ratio R. However, also note that a flat pion wave function was concluded in the framework 
of covariant quark-pion model j2J] and of QCD sum rules (2^]2(|. 

The above results seem not to be well consistent with the universality of the pion wave function. We emphasize that 
the infrared structures of the processes irj* — > 7 and n^/* — > ir are different at next-to-leading twist. For example, 
the collinear divergences associated with the pseudo-scalar structure 75 are absent in irj* — > 7, but exist in irj* — > ir. 
The three-parton wave functions contribute to both the pion transition form factor and the pion form factor. It is 
expected that at the maximal energy scales around 8 GeV 2 , where data are available, these subleading contributions 
are sizeable relative to the leading ones. It has been explicitly demonstrated that if higher-twist controbutions from 
parton transverse momenta are taken into account, the CZ wave function is not excluded by the data of the pion 
transition form factor |2^| , contrary to the conclusion in jis|| . It has been shown that the pion form factor suffers 
substantial higher-twist contributions proportional to the chiral condensate at currently available energy scales p8[ . 
The B — > ir transition form factor also receives nonvanishing higher-twsit contributions |2^] . Because of these next-to- 
leading-twist ambiguity, we argue that the above different extractions of the leading-twist pion wave function should 
not be regarded as an inconsistency. 
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Figure Captions 



1. Fig. 1: Lowest-order diagrams for 7r7* — > 7 (B — > 7^), where the symbol x represents the virtual photon (weak 
decay) vertex. 

2. Fig. 2: 0(a s ) radiative corrections to Fig. 1(a). 

3. Fig. 3: (a)-(d) Infrared divergent diagrams factored out of Fig. 2(d) and 2(e). (e) The graphic definition of the 
leading-twist pion wave function. 

4. Fig. 4: (a) Factorization of 0(a^ +1 ) reducible diagrams, (b) The Ward identity, (c) Factorization of 0(a^ +1 ) 
irreducible diagrams. 

5. Fig. 5: Lowest-order diagrams for ttj* — > n (B — > irW), where the symbol x represents the virtual photon 
(weak decay) vertex. 

6. Fig. 6: 0(a s ) radiative corrections to Fig. 5(a). 

7. Fig. 7: 0(a s ) radiative corrections to Fig. 5(b). 

8. Fig. 8: (a)-(d) Infrared divergent diagrams factored out of Fig. 2(d) and 2(e) for the B — > decay, (e) The 
graphic definition of the leading-twist B meson wave function. 
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